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ABSTRACT 

Context. Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By 
combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars 
in a few hours of observations. 

Aims. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, 
and estimate the dynamic range from interferometric observations. 

Methods. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools 
that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit using all inter¬ 
ferometric observables, i.e., the squared visibilities, closure phases and bispectrum amplitudes. The search procedure is made on a 
N X N grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the 
companion in the number of sigmas. The code CANDID also incorporates a robust method to set a 3cr detection limit on the flux 
ratio, which is based on an analytical injection of a fake companion at each point in the grid. Our injection method also allows us to 
analytically remove a detected component to 1) search for a second companion; and 2) set an unbiased detection limit. 

Results. We used CANDID to search for the companions around the binary Cepheids VI334 Cyg, AX Cir, RT Aur, AW Per, SU 
Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting VI334 Cyg and AX Cir were detected 
at > 25cr and > 13cr, respectively. The astrometric positions and flux ratios provided by CANDID for these two stars are in good 
agreement with our previously published values. The companion around AW Per is detected at more than 15(T with a flux ratio of 
/ = 1.22 ± 0.30 %, and it is located at p = 32.16 ± 0.29 mas and PA = 67.1 ± 0.3°. We made a possible detection of the companion 
orbiting RT Aur with / = 0.22 ± 0.11 %, and at p = 2.10 ± 0.23 mas and PA = -136 ± 6°. It was detected at 3.8cr using the closure 
phases only, and so more observations are needed to confirm the detection. No companions were detected around SU Cas and T Vul. 
We also set the detection limit for possible undetected companions around these stars. We found that there is no companion with 
a spectral type earlier than B7V, A5V, FOV, B9V, AOV, and B9V orbiting the Cepheids VI334 Cyg, AX Cir, RT Aur, AW Per, SU 
Cas, and T Vul, respectively. This work also demonstrates the capabilities of the MIRC and PIONIER instruments, which can reach a 
dynamic range of 1:200, depending on the angular distance of the companion and the (u, v) plane coverage. In the future, we plan to 
work on improving the sensitivity limits for realistic data through better handling of the correlations. 
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1. Introduction 


Send offprint requests to'. A. Gallenne 


Long-baseline interferometry (LBI) enables us to spatially re¬ 
solve components in close orbits, providing astrometric positions 
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at < 50 milli-arcsecond (mas) scale with micro-arcsecond accu¬ 
racy. When combined with spectroscopic radial velocities, we 
can obtain model independent estimates of the stellar masses and 
orbital parallaxes, which are fundamental parameters that help 
us study stellar properties and evolution. However, so far LBI 
is limited to bright stars (H < 1 mag) with typical magnitude 
differences AH < 6 mag. Adaptive optics imaging with single¬ 
dish telescopes reaches better contrasts, down to AH ~ 12 mag 
( Zurlo et al.|2014| l, but the angular resolution is limited to 0.2" at 
these detection levels. Long-baseline interferometry is therefore 
a complementary technique in terms of spatial scale by probing 
the innermost regions. 


LBI can detect components down to ~ 1 mas in the infrared, 
but the main limitation is sensitivity to high-contrast compan¬ 
ions. We roughly know the performances reachable by current 
interferometric combiners. So far, the faintest companion de¬ 
tected with LBI has a flux ratio of 0.75 % and orbits a Cepheid 
star of magnitude H - 3.85 mag ( Gallenne et al.|20T^ . Objects, 
such as brown dwarfs and hot giant planets, are still inaccessible 
because of a lack of sensitivity and accuracy of the instruments. 
Absil et al. (2011|l demonstrated a possible dynamic range of 


1;500 with the VLTI/PIONIER instrument (Le Bouquin et al. 
|2011| l, but this range has not yet been reached. To achieve this 
detection level, several hours of multitelescope observations are 
required to obtain as many simultaneous interferometric mea¬ 
surements as possible. 


Deriving the sensitivity limit from imaging can be deter¬ 
mined directly from the noise level, however, this is not the 
case for interferometric observations. There are some papers in 
the literature that discuss detection limits and methods to search 
for companions from interferometric data 
Bouquin & Absil|2012| l, but these studies 
mgs: the searching method is not formalized, the sigma detec¬ 
tion is not robust, and they do not take the bandwidth smear¬ 
ing into account. Therefore, a robust implementation to search 
for components does not exist thus far. This kind of method is 
particularly critical to detect faint companions, as they can be 
at the sensitivity limit of the instrument or even not be a sta¬ 
tistically signihcant detection (i.e., < 3cr). We therefore created 
CANDID (Companion Analysis and Non-Detection in Interfero¬ 
metric Data) to address these aspects. This is a suite of Python 
tools, which contains two main functions: 1) one to perform a 
systematic search for faint companions (Sect. [^; and 2) one to 
estimate the detection limit from long-baseline interferometric 
observations (Sect. 0. This tool is made available to the com- 
munitjQ 

CANDID is made for detecting high-contrast, point-source 
companions orbiting a spatially resolved primary star, although 
it also works with an unresolved primary and contrast <50%. 
In this paper, we used CANDID to look for companions in binary 
Cepheids. We present the first main function to search for com¬ 
panions in Sect. verify and clarify the detection level of our 
previously detected faint companions, and we report new detec¬ 
tions for other Cepheids. In Sect. we present the second main 
function of CANDID and explain our robust method to set detec¬ 
tion limits from interferometric data. In Sect. we then use this 
formalism to our set of Cepheid observations to derive the detec¬ 
tion limits. We finally present our conclusions in Sect.j^ 


( Absil et al.|201 1 Le 


have some shortcom- 



2. Searching for companions 

Our grid search is traditionally performed in a three-dimensional 
space. Specihcally, we vary the position of the component 
{Aa, A6) and the companion/star flux ratio /, and then compute 
the for each of these positions and flux ratios (the third di¬ 
mension being /). The weakness of this method is the resolution 
of the grid, i.e., if the grid is to too coarse, the detection can be 
missed. 

In addition, for high-contrast binaries and/or low accuracy 
data, s.x^ map can show fake or nonsignificant detections. Some¬ 
times, some authors show the estimated x^ map with the most 
probable location of the component without any detection level 
mentioned, only the reduced x^ variation is given. This parame¬ 
ter, however, is not optimal for checking the detection level be¬ 
cause it depends on the number of degrees of freedom (dof). The 
large quantity of data required to detect a faint companion using 
LBI results in a large dof. For instance, we can have the most 
probable location with the lowest reduced chi-square, x^, equal 
to 1.0, and the highest in the whole map equal to 1.1, i.e., with 
only a variation of 0.1. With this information, however, there is 
no way to know if the detection is statistically significant or not. 
In this section, we therefore present a systematic approach to 
search for components and set the detection level using the x^ 
and the number of dof. We then applied it to the case of binary 
Cepheids. We verify our previous detections for the Cepheids 
V1334 Cyg and AX Cir hrst ( [Gallenne et al.|2013b|[20T4) l, and 
then report new detections for RT Aur and AW Per. 


2.1. Detection method 


A more rigorous approach is to perform a grid of fit using a least- 
squares algorithm, with a starting grid spacing that is guaranteed 
to find the global minimum. The grid in question is the 2D grid 
of starting points for the companion position. For each starting 
position, a multiparameter fit is performed: the companion posi¬ 
tion and its flux ratio (possibly the stellar diameters) are adjusted. 
Each position of the grid leads to a local minimum. Ideally, if the 
starting grid is fine enough, multiple starting points lead to the 
same local minima: this guarantees that all the local minima are 
explored and that the global minimum is indeed the global min¬ 
imum. Hence, the criteria to decide if the global minimum is 
global is a posteriori. In CANDID, we require (a posteriori) that 
on average, each unique minimum is reached from two starting 
points of the starting grid. We also provide statistics on the "trav¬ 
eling" distance of the fit (between the starting and end points), 
compared to the size of the grid. We require that the median 


travel distance should be less than V2/2 of the size of the square 
grid. The 10 and 90 % percentiles are also provided to the user 
to assess the typical travel distance. 

A systematic search using a grid of fit is an iterative process. 
First, a coarse starting grid is chosen and the fits are performed. 
CANDID estimates, based on the traveling distances of the hts, 
how to refine the starting grid. A second series of fits are run, 
using the finer grid, which might take much longer, but the global 
minimum can be trusted to be global (at least within the area 
searched). 

We searched for components with a maximum distance to 
the main star of 50 mas. For a wider range, the loss of coherence 
caused by spectral smearing of the companion is the main limita¬ 
tion because it degrades the dynamic range. A complete discus¬ 
sion about the search region is presented by [Absil et ah] ( 201 l| l 
and [Le Bouquin & Absil[ ( |2()12[ ). The main limitation is the spec¬ 
tral sampling compared to the relative position of the component. 
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i.e., to avoid a significant smearing we need R > pB^jA, with R, 
Bp and p the spectral resolution, the projected baseline, and the 
separation, respectively. For PIONIER, with R = 18 and a mean 
projected baseline Bp = 100 m (used for our observations), we 
have p < 60 mas, and for MIRC with R - 42 and B = 200 m, we 
have p < 70 mas. [Zhao et al!] ( |2007| l recommended a more strin¬ 
gent criteria p < RAjiSB) to assure uncorrupted data, leading 
to p < 15 mas for PIONIER and MIRC. However, a companion 
can still be detected at separations larger than 15 mas by tak¬ 
ing bandwidth smearing effects into account. We chose 50 mas 
as our limiting range to avoid making the grid search too strin¬ 
gent. In addition, companions located at more than 50 mas are 
detected more efficiently using adaptive optics on a single-dish 
telescope (through imaging or sparce aperture masking). 

Each point in the grid is fitted with the following binary 
model, representing a spatially resolved primary star with a point 
source component: 


V{u, v) 


Vi,{u,v) + G{OfVc{u,v) 

1 +/ 


with. 


Vi,iu,v) 


2J\ (x) 

X 


( 1 ) 


( 2 ) 


Vc(u, v) = exp[-2i7r(MAa -H vA(5)//l], 


( 3 ) 


G(0- 


sinf 




with ^ = 


niuAa + vA5) 
Rd ’ 


( 4 ) 


where 7i(x) is the first-order Bessel function, x = 
n6\jo -H v^/d, (m, v) the spatial frequencies, 0 ud the uniform 
disk angular diameter of the primary star, d the wavelength, 
/ the flux ratio between the companion and the primary star, 
(Aa, AS) the relative position of the component with respect to 
the primary, and R - AjAA the spectral resolution. The function 
G(^) is a corrective term to overcome the effect of bandwidth 
smearing ( Lachaume & Berger|2013| l. The fitted parameters are 
0UD, /, Aa, and AS. Although the bandwidth smearing can be fit 
in CANDID, we kept it fixed for this work because we noticed 
that it cannot be constrained by these data. 

The interferometric observables, i.e., the squared visibility 
y^, the closure phase CP, and the bispectrum amplitude Bamp, 
are then estimated from the squared modulus and the bispectrum 
in closed triangles as follows: 


y2 = |y|2 and 


( 5 ) 


which provides the bispectrum amplitude and the closure phase 
from the definition B - 


Bamp = |B| and CP = arg(B). (6) 

The data are then fitted simultaneously using a Levenberg- 
Marquardt least-squares minimization algorithm with 

^ (CBo - CBni)Vo-cp + ^ (Bamp.o “ Bainp,m)^/o-B^p 

+ ( 7 ) 

where the indexes o and m denote the data and the model, respec¬ 
tively. CANDID can fit all observables or just one, depending on 
the data available and the user. We then divided by the number 


of degrees of freedom to obtain a map of the chi-square minima, 
and interpolated these minima on a regular grid. 

The grid resolution is critical as it depends on the distance 
explored from the initial and final positions. This is why we 
implemented an estimate of the optimum grid resolution in our 
code. It is worth mentioning that a denser starting grid is not nec¬ 
essarily better as the fit would not be improved, and we would 
lose in computation time (CANDID performs a 30 x 30 grid with 
1 mas steps in 50 s using six cores and a data set with 879 de¬ 
grees of freedom). To make the grid search faster, CANDID was 
developed for parallel processing on multicore machines. 

Once the map of the minima is computed, we can check the 
variation in the whole map, however this does not tell us about 
the statistical significance: specifically, whether the most prob¬ 
able location is detected at Icr or more. Although it might not 
be important for low-contrast companions for which the varia¬ 
tion is large enough, this is critical for components with a flux 
ratio < 5 % with a variation of the minima of a few percent. The 
number of dof is an important parameter in that context. Assum¬ 
ing that the data follow Gaussian statistics, we implemented in 
CANDID an estimate of the number of sigma for each point in the 
grid in order to obtain a ncr detection map. The formalism we 
used is based on the probability P (or confidence interval) with v 
degrees of freedom, as already employed by Absil et al. ( 2011) 1. 
The number of sigma demonstrates how our binary model is sig¬ 
nificant compared to a uniform disk model (i.e., a single star). 
We used the following formula for the probability: 


P{Aa, A6)^\- CDFy 


kTuD 




( 8 ) 


where and are the minimum chi-square for the binary 
and the uniform disk model (i.e., fitting eq. [T] and [^ respec¬ 
tively), and CDF denotes the cumulative probability distri¬ 
bution function with v degrees of freedom. We then convert the 
probability into the number of sigmas, ncr (e.g., 99.73 % = 3cr, 
99.99 % = 4cr, ... ). To avoid big float numbers, we limited the 
maximum value to 50cr. 

This formalism therefore provides map to find the most 
probable location of a companion, if any, and a ncr map giving 
the detection level at each point in the grid. 


2.2. Previously published detections 

To validate our method, we computed the maps of the x^ min¬ 
ima and estimated the detection level for our previously de¬ 
tected companions around the Cepheids VI334 Cyg and AX Cir 
([Gallenne et al.||2013b| [2014) observed with the instruments 
CHARA/MIRC and VLTI/PIONIER, respectively). The journal 
of these previous observations are repotted in Table)!] The maps 
are presented in Eig. □ and [^ for which all observables were 
fitted, except for PIONIER where only CP and were used be¬ 
cause there is no good estimator of B^mp so far. The central part 
has been hidden to improve the clarity of the intensity map level, 
which can be biased by the primary star. The companion orbit¬ 
ing V1334 Cyg is detected at more than 25cr, and at more than 
13cr for AX Cir. We summarized the detection levels for these 
two stars in Table [^ including fitting only the closure phase sig¬ 
nal. The CP is more sensitive to faint off-axis companions and 
is also less affected by instrumental and atmospheric perturba¬ 
tions than the other observables (i.e., and B^mp). Pitting all of 
the observables can improve the detection level because we add 
more information, but it can also affect the results, depending on 
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Table 1. Journal of the observations. 


UT 

MJD 

Star 

^bracket 

Configuration 

Nv2 

Ncp 

Calibrators 

2012 Jul. 26 

56135.449 

V1334Cyg 

2 

S1-S2-E1-E2-W2 

48 

42 

HD 200577, HD 214200 

2012 Sep. 30 

56201.221 

VI334Cyg 

3 

S1-S2-E1-E2-W1-W2 

62 

68 

HD 185395, HD 199956, 
HD 218470, HD 207978 

2013 Jul. 14 

56487.983 

AX Cir 

6 

DO-HO-Gl-Il 

300 

200 

HD 133869, HD 129462 

2012 Jul. 26 

56134.354 

T Vul 

2 

S1-S2-E1-E2-W1-W2 

199 

260 

HD 192518, HD 205852 

2012 Sept. 30 

56200.228 

T Vul 

2 

S1-S2-E1-E2-W1-W2 

175 

210 

HD 189849, HD 198692, 
HD 207978 

2012 Sept. 30 

56200.434 

SU Cas 

4 

S1-S2-E1-E2-W1 

259 

255 

HD 12216, HD 19267, 
HD 34200 

2012 Oct. 01 

56201.428 

AW Per 

1 

S1-S2-E1-E2-W1-W2 

210 

280 

HD 19845, HD 30825, 
HD 35940 

2012 Oct. 01 

56201.507 

RT Aur 

1 

S1-S2-E1-E2-W1-W2 

105 

140 

HD 48682 


Notes. Nbrackct: number of data blocks. Nyi and Ncp: number of squared visibilities and closure phase. Adopted calibrator diameters: HD 200577 
= 0.758 ± 0.052 mas, HD 214200 = 0.790 ± 0.050 mas, HD 185395 = 0.726 ± 0.014 mas, HD 199956 = 0.603 ± 0.043 mas, HD 218470 = 0.477 ± 
0.033mas, HD 207978 = 0.571 ±0.040 mas, HD 133869 = 1.043 ± 0.015 mas, HD 129462 = 0.857 ±0.061 mas, HD 192518 = 0.418 ± 0.029 mas, 
HD 205852 = 0.461 ±0.032mas, HD 189849 = 0.510±0.036mas, HD 198692 = 0.660±0.056mas, HD 207978 = 0.571 ±0.040mas, HD 12216 = 
0.467±0.033 mas, HD 19267 = 0.586±0.042 mas, HD 34200 = 0.652±0.046mas, HD 48682 = 0.616±0.043 mas, HD 19845 = 0.788±0.056mas, 
HD 30825 = 0.564 ± 0.040 mas, HD 35940 = 0.615 ± 0.044 mas 



Aa (mas) 
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Fig. 1. map of the local minima (left) and detection level map (right) of V1334 Cyg for the observations on 2012-10-01. The yellow lines 
represent the convergence from the starting points to the final fitted position. The maps were reinterpolated in a regular grid for clarity. The axis 
limit was chosen according to the location of the companion. 


the magnitude of the biases. We notice from Table|^that includ¬ 
ing and Bamp degrades the detection level for VI334 Cyg, 
although it is still significant, and including these observables 
improves the detection for AX Cir. A possible explanation that 
addresses why the detection level decreases sometimes when we 
add the is that we add correlated noise, which is not consistent 
with our hypothesis of uncorrelated noise. When only the CP is 
used, the angular diameter of the primary is first determined by 
fitting only a uniform disk model to the square visibilities, and 
then kept fixed during the grid search. 


The resulting fitted parameters (i.e., the astrometric position, 
flux ratio, and the angular diameter) are in good agreement with 
the values determined in |Gallenne et al.| ( |2()13b|[2014| l. 
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2.3. New detections of the companions around the Cepheids 
RT Aur and AW Per 


Our interferometric program on Galactic binary Cepheids, 
which started two years ago, is promising to directly mea¬ 
sure the dynamical masses. We have obtained several observing 
nights in both hemispheres with the multitelescope combiners 
CHARA/MIRC and VLTI/PIONIER to detect the close compan¬ 
ions of a few Cepheids ( [Gallenne et al.|20T4 |2013b|a l. Here we 
report new detections for the Cepheids AW Per and RT Aur. 


The observations were performed in 2012 using the Michi¬ 
gan InfraRed Combiner (MIRC) installed at the CHARA Array 
( |ten Brummelaar et al.|2005l l, located on Mount Wilson, Califor¬ 
nia. The CHARA Array consists of six 1 m aperture telescopes in 
a Y-shaped configuration (two telescopes on each branch), ori¬ 
ented to the east (El, E2), west (W1,W2) and south (SI, S2), 
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Fig. 2. Same as Fig.[^ except for AX Cir for the observations on 2013-07-14. 


providing good coverage of the (u, v) plane. The baselines range 
from 34 m to 331 m, providing high angular resolution down to 
0.5 mas at H. The MIRC instrument ( |Monnier et al.|2004||2010| l 
is an image-plane combiner, which enables us to combine the 
light coming from all six telescopes in H or K. MIRC also offers 
three spectral resolutions {R - 42,150 and 400), which provide 
15 visibility and 20 closure phase measurements across a range 
of spectral channels. 

Table 2. Summary of the detection level for our previous published data 
and new detections. 


Star 

CP only 
(no-) 

ALT'* 

(no-) 

dof 

Instrument 

VI334Cyg 

42.1 

25.4 

1583 

MIRC 

AX Cir 

13.6 

19.8 

1499 

PIONIER 

RT Aur 

3.8 

2.3 

879 

MIRC 

AW Per 

19.6 

15.0 

1757 

MIRC 


Notes. All means CP +V'^ + Samp, and Sa^p is only used with MIRC 
data. 


We observed the short-period Cepheids RT Aur (HD 45412, 
Ppuk = 3.73 d) and AW Per (HD 30282, = 6.46 d) with all 

six telescopes. We used the //-band filter with the lowest spec¬ 
tral resolution in which the light is split into eight spectral chan¬ 
nels. Table lists the journal of our observations. We followed 
a standard observing procedure in which we monitored the in¬ 
terferometric transfer function by observing a calibrator before 
and/or after the Cepheids. The calibrators, listed in TablefTl were 

_ . — m ~ _ J— f 


select ed using the 5earc/iCa0 software (Bonneau_ et al. 


201l| provided by the Jean-Marie Mariotti Centeij^ 


2006 


We reduced the data using the standard MIRC pipeline 
( [Monnier et al.|2007| ), which consists of computing the squared 
visibilities and triple products for each baseline and spectral 
channel, and correcting for photon and readout noise. Squared 
visibilities are estimated using Fourier transforms, while the 
triple products are evaluated from the amplitudes and phases be¬ 
tween three baselines forming a closed triangle. 


^ Available at http://www.jmmc.fr/searchcal 
^ http://www.jmmc.fr 


We used CANDID to search for a component within +50 mas. 
For RT Aur, we might have detected a companion at 3.8cr using 
only the closure phases, while we only have a 2.3cr detection 
using all observables. The possible companion is detected at p = 
2.1 mas and PA - -136° with a flux ratio / = 0.22%. The 
component orbiting AW Per is detected at > 15cr at p = 32 mas 
and PA - 67° with a flux ratio / = 1.22%. The grids of fit 
and detection level maps are shown in Fig and and the final 
fitted parameters are listed in Tabl e For the u ncertainties, we 
used the conservative formalism of [Boffin et al.| ( |20T4j l for all the 
fitted parameters, i.e., 

(t\ = //spO-L + O.OOOIA^ (9) 

where Ajp is the number of spectral channels and X denotes the 
fitted parameters (i.e., Ax,Ay,f and 0 ud)- The first term takes 
into account that the spectral channels are almost perfectly cor¬ 
related, and the second term comes from the fact that the wave¬ 
length calibration is only precise at a 1 % level. The parame¬ 
ter (Tstat is the statistical error from the bootstrapping technique 
(bootstrap on the calibrated data with replacement) using 10 000 
bootstrap samples (also included in CANDID). We then took from 
the distributions the median and the maximum value between 
the 16% and 84% percentiles as the uncertainty (although the 
distributions were roughly symmetrical). 

There are additional significant peaks in the RT Aur maps, 
with a detection level > 3cr, however, these are spuriously pro¬ 
duced by the (m, v) coverage and the presence of the companion. 
It is worth mentioning that RT Aur was only observed for one 
hour (one sequence), and we need more data to confirm the pres¬ 
ence of the companion. A more complete discussion about the 
detected companions is presented in Sect. 

2.4. Undetected companions 

We also observed with MIRC the Cepheids SU Cas (HD 17463, 
P = 1.95 day) and T Vul (HD 198726, P = 4.44 day). We used 
five and six telescope configurations, using the same instrument 
setup and calibration procedures as explained in the previous 
section. Table[^lists the journal of these observations. 

We did not detect any significant companions around these 
stars, either using all of the observables or only the CP, i.e.. 
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Fig. 3. Same as Fig.[^ except for RT Aur using only the closure phase. 

lo.go 

^ 0.96 

- 0.93 
= 0.90 'i' 

m 0.87 ^ 

- 0.84 

- 0.81 
- 0.78 

40 30 20 10 0 -10 -20 -30 -40 

An (mas) 

Fig. 4. Same as Fig.[T| except for AW Per. 

Table 3. Final hest-fit parameters. 



40 30 20 10 0 -10 -20 -30 -40 


An (mas) 




RT Aur 
(p = 0.32 

AW Per 
(p = 0.52 

SU Cas 
(P = 0.77 

T Vul 

(p ^ 0.27 1 (p = 0.12 

6»ud (mas) 

0.699 + 0.011 

0.627 + 0.018 

0.609 + 0.043 

0.608 + 0.013/0.635 +0.018 

/(%) 

0.22 + 0.11 

1.22 + 0.30 

- 

- 

Aa (mas) 

-1.458 + 0.238 

29.624 + 0.305 

_ 

- 

Ad (mas) 

-1.506 + 0.224 

12.523 + 0.147 

- 

- 


Notes. 0UD: uniform disk angular diameter, respectively. /, Ax, Ay: flux ratio and position of the companion. 


no more than 1.9cr and 2.9cr for SU Cas and T Vul, respec¬ 
tively. However, we were able to measure their angular diam¬ 
eters, which are listed in Table A more detailed discussion of 
these binary systems is presented in Sect. 


3. Detection limit of high-contrast binaries 


It is important to be able to check the dynamic range that can 
be reached with a given set of data and any interferometric com¬ 
biner. Absil et al. ( |201 l| l already presented a method to set de¬ 
tection limits for the VLTI/PIONIER instrument, but it has some 
shortcomings, that is why we propose a more robust formalism. 


Article number, page 6 of 


12 




























































A. Gallenne et al.: Robust high-contrast companion detection from interferometric observations 



Fig. 5. Same as Fig|^ except we analytically removed the companion. 


Absil’s method: Their method is based on comparing a uniform 
disk model with a binary model for each position (Aar, A6) in the 
grid. They then checked whether the probability of the binary 
model is consistent with the data using. 


where we kept the term to avoid having an (small) offset (see 
next section). For a high-contrast companion, i.e., for / << 1, 
we can approximate the bispectram at the first order in / as. 


P( Aa, A5) = 1 - CDFy 


^ Tud j 


( 10 ) 


We notice that this equation has a different ratio in the CDF than 
Eq. This is because we assumed that the binary model is the 
trae model, while |Absil et al.| ( |2011| l assumed the uniform disk 
as the true model. In theory, both equations should lead to the 
same results, however, as we see in subsequent sections, their 
method is more sensitive to biased data and can sometimes lead 
to under- or overestimated detection limits. 


Our method: We suggest an alternative method, which is based 
on the injection of a companion into the data at each astromet¬ 
ric position with different flux ratios. As we inject a compan¬ 
ion, we therefore know that the binary model should be the true 
model, and we can use Eq. |^to obtain the probability of the bi¬ 
nary model to be the trae model. We introduced this method be¬ 
cause we think it is more robust, as we demonstrate in the next 
section. 

In this section, we first introduce approximate formulae for 
high-contrast companions. We then explain how we inject an ad¬ 
ditional component into the data and derive the detection limits. 

3.1. High-contrast approximation 

The complex visibility for a binary system composed of a re¬ 
solved primary star and an unresolved component is given by 
Eq. [T] The squared visibility is given by. 


1 -H/^ 

^ Vl+2G(OfV^ COS ip+ G(0^f 

( l +/)2 


( 11 ) 


B = 


(Vi + G(A)/e‘‘^‘)(V2 + G{C2)fe'‘f^-){Vi + G((2)fe-^‘P^) 
( 1 +/)' 




( 1 +/)^ 


1 +/ 


Eh 


E2* 


V: 


3* 


.z 


( 12 ) 


where B* = 

The bispectrum amplitude and the closure phase are then es¬ 
timated as, 

Bamp = ^y^^|B*||Z| (13) 

CP = arg(B*) - arg(Z). (14) 


3.2. Adding/removing a component 


If no companion is detected in the data, we can assume that the 
measured values only represent the primary star, i.e., a uniform 
disk, and we can substitute the index in the previous equa¬ 
tions by the index "obs" (note that the function Z also depends 
on the primary star visibility). It is now simple, from Eq. 
and[^ to inject a companion to the observed data, which 
sponds to the following equations for the following new observ¬ 
ables: 


TT][T3 

corre- 


y2 

+ 2G{OfVobs cos ip + GiO^f 


(l+/)2 


^amp 

“ + _/)^ ^amp,obs I'^li 

(16) 

CP 

= CBobs - arg(Z). 

(17) 


The original "oifits" files provide all the necessary spatial 
and spectral information to reconstruct the individual phases and 
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visibilities. We also take the coherence loss effect caused by the 
spectral smearing of the companion into account. 

Inversely, if a component is detected, we can use this for¬ 
mula to analytically remove the companion (with a negative flux 
ratio), and then check for another possible fainter component or 
estimate the detection limit to rule out any other companions. 
This step is critical to obtain unbiased detection limits. As an 
example in Fig. , we analytically removed the detected com¬ 
ponent orbiting RT Aur, and notice that there is no significant 
detection of a third component and that the other minima are 
not additional companions. 

To check our approximation, we created a model of a single 
star with a uniform disk of 1 mas (3h of PIONIER observations 
with a point every 30 min), for which we added a companion at 
a position Aa = Ab = 50 mas with a flux ratio f - 5%. We then 
compare this model with a true binary model using the same pa¬ 
rameters (Eq.0. The difference between the approximation and 
the true model is shown in Eig. We found for the amplitude of 
the bispectrum a standard deviation of the relative error < 0.3 %, 
and a difference < 0.005° for the closure phase, which is lower 
than the achievable interferometric accuracy. Closer companions 
with the same or higher contrast (i.e., < 5%) lead to smaller er¬ 
rors. 

We kept the term in the squared visibility, otherwise it 
gives an offset in the relative error of about 0.3 % for f - 5%, 
and decreases with decreasing flux ratios. Although this value 
is negligible compared to the current possible data accuracy 
(AV^IV^ = ABamp/Bamp ~ 2%, ACP ~ 0.5°), the offset is 
larger for lower contrasts (~ 1.2% with / = 10%). As no ap¬ 
proximation is made for V^, the formula is valid for any flux 
ratio. We therefore also checked the bispectrum approximation 
for brighter components, and found that we can use it up to a 
flux ratio of 50 % only using CP and V^. We listed in Table 
the relative error for different flux ratios and various positions 
(1 < p < 50 mas and 0 < PA < 2n). We see that CANDID can 
also be used for low-contrast companions. 



B/A (1/rad xlO^ ) 

Fig. 6. Difference between the true binary model and the approximation 
model. For the bispectrum, the abscissa is the maximum of the three 
spatial frequencies. 


3.3. Estimating the detection iimit 

We implemented in CANDID a tool to estimate the detection limit 
for a given set of interferometric data. The method is to compare 
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Table 4. Validity domain of our approximation. 


/ 

(%) 

ACP 

n 

^^amp /^amp 
(%) 

5 

<0.01 

< 0.5 

10 

<0.02 

< 1.9 

50 

<0.25 

< 67 

75 

< 0.7 

> 100 

100 

< 3.5 

> 100 


the x^ obtained for a model without a companion and the x^ ob¬ 
tained for a model with an injected companion. The method does 
not work, a priori, with an already detected companion, first be¬ 
cause the previous approximation might not be valid and also 
because the resulting detection limit would be biased by a sys¬ 
tematic offset related to the flux ratio between the components. 
This means that any detected component has to be analytically 
removed first. 

The method is based on an A x A grid with a range of 
+50 mas, for which the minimum A depends on the optimum 
resolution estimated from the^^ map (see Sect.|^. At each point 
in the grid, we inject a companion with various flux ratios and 
we compute thex^- As we know that the true model is the binary 
model (because we injected a companion), we used our previous 
equation, Eq. to estimate the number of sigmas for each flux 
ratio. We then interpolated the flux ratio values at 3cr, which we 
set as the significance level. This means that lower flux ratios are 
not detected significantly. Doing this for all points in the grid, 
we then have a 3cr detection limit map for the flux ratio. To have 
a quantitative estimate of the sensitivity limit with respect to the 
separation, we estimated a radial profile, / 3 o-(r), using the 90 % 
completeness level (i.e., 90 % of all possible positions) from the 
cumulated histogram in rings for all azimuths. This tool also in¬ 
cludes parallel processing to make the calculation faster. 

In theory, for uncorrelated data with Gaussian statistics, our 
method should lead to the same results as |Absil et al.| ( |20TT| l in 
terms of the detection limit. However, real data are often biased 
by different sources (atmospheric turbulence, mechanical vibra¬ 
tions, ...). We performed two tests to compare both methods, one 
including uncorrelated Gaussian noise and another with a noise 
model, i.e., correlated non-Gaussian noise. We used all observ¬ 
ables and three data sets for each test. The Aspro2 softwar^ 
were used to create these synthetic data sets, however, it does not 
have implemented the bandwidth smearing effect. We therefore 
did not take it into account in CANDID (this would not change the 
conclusion of our test). 

First test: We created the first "ideal" synthetic data set (i.e, 
without noise) representing a uniform disk of 1 mas (3h observa¬ 
tion with three spectral channels with the PIONIER instrument). 
We then added uncorrelated Gaussian noise and estimated the 
detection limit using the Absil’s method and our formalism. As 
expected, we see in Fig. that both methods provide the same 
results. We then created a second data set by adding a component 
at Aa - A6 - 5 mas with / = 2 % as a source of bias to the pre¬ 
vious model. The companion creates a departure from the ideal 
measurements (like noise) but in such a way that the observa¬ 
tions at different spatial frequencies have a correlated departure 
(hence, correlated noise). We then estimated the detection lim¬ 
its for this second data set, and although the trend between the 
two methods are similar, there is an offset of ~ 2 % because of 

Available at http: / /www . jimnc. fr/aspro_page. htm 
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Fig. 7. Comparison between the CANDID detection limit tool and the 
formalism of |Absil et al.H20lT| l for a uniform disk model with Gaussian 
noise. 


the presence of the faint companion. Finally, we analytically re¬ 
moved the companion to get our third data set. Comparing the 
three data sets, i.e., estimating the total variation of the detec¬ 
tion limit for each radius, we noticed that both methods vary in 
a similar way. This was expected as all of the noise sources in¬ 
serted have Gaussian statistics. However, it is worth mentioning 
that the scatter tends to be a bit larger for the Absil’s method for 
increasing flux ratios and wider separations of the companion 
(due to the correlated noise introduced by the component). 


Second test: We created two new data sets similar to the first 
test, i.e., one for a uniform disk of 1 mas and another for a bi¬ 
nary system with Aa = A6 - 5 mas and f - 2% (with average 
atmospheric conditions, which is an option for the noise model). 
The only difference is that the noise is no longer Gaussian; it is 
represented by a more complicated model that includes the in¬ 
strument response, atmospheric turbulence, photon, and detector 
nois^ As shown in Fig. the absil’s method is more sensi¬ 
tive to the presence of the companion, which in this case gives 
a higher sensitivity limit. The third data set is created by remov¬ 
ing the component using our approximation. We then compared 
the variation of the detection limit between all of the data sets, 
as described previously, i.e., the minimum and maximum value 
between the three synthetic data sets at each separation. Fig. 
shows that the scatter from the Absil’s method is larger. We per¬ 
formed additional tests with lower flux ratios and at different 
positions, and the scatter from the Absil’s method is larger than 
our formalism most of the time. We therefore conclude that our 
method is more robust to biased data when estimating detection 
limits. 

In our tests, the limits derived from the injection method are 
lower than when we use Absil’s method, but this is not always 
the case for real data. The injection method occasionally results 
in higher sensitivity limits, depending on the magnitude of the 
biases. We find that Absil’s method may under- or overestimate 
the detection limits depending on the data set. However, in some 
cases when the data are not affected by significant biases, both 
methods give similar results. 


^ Details are explained in http://www.jmmc.fr/doc/approved/ 
JMMC-MEM-28S0-SO01.pdf 



Fig. 8. Comparison between the CANDID detection limit tool and the 
formalism of |Absil et al.H20fT} for a uniform disk model with a noise 
model, biased by the presence of a faint companion. 



radial distance (mas) 

Fig. 9. Total variation (minimum and maximum values) between the 
three synthetic data sets with a noise model for both methods: a single 
star, a single star -t- an unresolved component as source of bias, and the 
previous model with the component analytically removed. The variation 
has been normalized to unity. 



Fig. 10. Flux ratio detection limit at 3cr for a second companion around 
V1334 Cyg (for 2012 Oct. 01). 

4. Detection limit of our binary Cepheids 

In this section, we set the detection limit for the sample 
of Cepheids previously presented. For VI334 Cyg, AX Cir, 
RT Aur, and AW Per, where a companion is detected, we re¬ 
moved analytically the companion first. Although both methods 
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rived from our injection method. We present two detection limits 

mu - 3.94 + 0.01 mag from Monson &. Kercel^O 

T]l. Using the 

for each star, one using all of the observables, and one using only 

/T-band P-L relation to get d = 428 pc (Storm et al. 

201 l|l, we es- 

the closure phases. We listed three different values, the average 

timate its spectral type to be later than an El star. T 

lis is compat- 


for r < 50 mas, r < 25 mas and r > 25 mas, which can be rel¬ 
evant when the limit increases with r. All of the final detection 
limits are listed in Table where the values are conservative as 
they correspond to the mean plus the standard deviation for the 
given radius range. 

From an evolutionary timescale point of view, most of the 
companions should be stars close to the main sequence. We 
therefore set upper limits for the spectral type of the compan¬ 
ion assuming it is on the main sequence, based on their //-band 
luminosities. 


VI334 Cyg: This binary system contains a visual and a spec¬ 
troscopic component. While the visual companion is > 150 mas, 
the close component was spatially resolved using interferometry 
by |Gallenne et al.| ( 2013b| l, and has a flux ratio of ~ 3.1 % in H 


ible with the A4V star upper limit determined by Evans ( 1992a|l 


(more details on this companion are presented in Gallenne et al. 
2013b|l. There is no evidence of a third component so far. Af 


ter removing the close companion analytically, we estimated the 
dynamic range as explained in the previous section. Between the 
two epochs, we reached a maximum average sensitivity limit of 
/ 3 o. = 1.54 %. The limit is lower for the first epoch because the 
atmospheric conditions were better than the second epoch. The 
contrast upper limit at 3cr is shown in Fig. with the average 
values for the two epochs in Table|^ This converts to magnitude 
difference of Aniu > -2.5 log/ 30 - = 4.5 mag. Using the distance 
d - 683 pc from the /T-band P-L relation for first overtone (FO) 
pulsators ( |Bono et al.|2002|l and the average magnitudes of the 
Cepheid H - 4.66 mag ( |Cutri et al.||2003| l, we can exclude the 
presence of additional companions with a spectral type earlier 
than B7V stars. 


AX Cir: The spectroscopic companion was first detected by 
|Lloyd Evans] (|1982|l, and was spatially resolved for the first 


time with FBI by Gallenne et al.| ( |2014[ with a more detailed 
discussion about this companion). The average detection lim¬ 
its are listed in Table We reached a dynamic range of / 30 - = 
0.36%, corresponding to Amjh > 6.1 mag if another compan¬ 
ion is present. With d - 500 pc (from the //-band P-L relation 
of Storm et al.|2011[ ) and H - 4.66 mag for the Cepheid ( |Cutri| 
|et al. 2003| l, we can rule out any other component with a spectral 
type earlier than A5V star. 


RT Aur: The binary nature of this short-period Cepheid is still 
uncertain, however we might have detected it for the first time. 
Some authors suggested the presence of an early-type compan¬ 
ion ( |Janot-Pacheco|[T97^ |Balona|ri977| l, while others did not 
see evidence of an additional component ( |Harris|[T981[ |Gieren| 
|1985| l. [Leonard & Turner] ( |1986| ) summarized various studies 
from that time about the possible companion orbiting RT Aur. 
Based on spectra from the International Ultraviolet Explorer 
satellite (lUE), [Evans ( 1992a| l did not report any detection and 
showed that any main-sequence secondary has to be cooler than 
an A4 star. Radial velocity measurements do not show any vari¬ 
ations from orbital motion. Recently, [Turner et al. ( 2007| l stud¬ 
ied long-term photometric light curves and reported a sinusoidal 
trend consistent with a light travel time effect in the binary 
system. The companion we might have resolved has a flux ra¬ 
tio of 0.21 + 0.12%, i.e., Amu - 6.7 + 0.6mag, leading to 


from the International Ultraviolet Explorer (lUE) spectra. How¬ 
ever, we need more observations to confirm the existence of the 
companion as we are at the sensitivity limit of this data set. 

Removing this possible companion from our interferometric 
data, we estimated the detection limit. We reached a minimum 
dynamic range of 0.47 %, corresponding to Amu >5.8 mag. This 
is a 90% completeness azimuthal value, which is a conserva¬ 
tive value, and the component at a given position with a slightly 
higher contrast might be detected (as in this case). This limit al¬ 
lows us to exclude the presence of an additional companion with 
a spectral type earlier than FOV. 

AW Per: This Cepheid is a spectroscopic binary with an or¬ 
bital period of ~ 40 yr. First discovered by [Miller & Preston 
( |1964| l, it took several years to derive the first orbit from radial 
velocity data ( [Evans[[T989| l. It is likely that this companion is 
itself a binary, as argued by Evans[ ( T989| l; Evans et ah] ( 2000 ) , 
because the magnitude difference between the Cepheid and its 
companion is not consistent with equal masses and predictions 
from evolutionary tracks. Unfortunately, we do not have enough 
angular resolution and sensitivity to detect this third companion. 
The properties of the brightest companion were studied based on 
lUE spectra by Evans (1994 1995| l, who found its spectral type 
to be B8.3V. This is in agreement with our detection with a flux 
ratio f - 1 .22 ± 0.30 %, i.e., a spectral type in the range B6-B9V 
(using d = 853 pc from a P-L relation and K - 4.63 mag for the 
Cepheid). We estimated for the companion mu = 9.6 ± 0.3 mag 
(using for the Cepheid at our given phase mu = 4.84 + 0.01 mag 
from [Monson & Pierce[[2011[ ). [Massa & Evans[ ( [2008| l also de¬ 
termined the angular separation of the component for another 
epoch, which allows us later, with more astrometric points from 
interferometry, to estimate all the orbital elements, including the 
inclination and the semi-major axis (see also [Gallenne et al.[ 

[2DT^ . 

We derived a maximum sensitivity limit at 3cr of 0.62 %, i.e., 
Aniu >5.5 mag. We can therefore exclude any additional com¬ 
ponents with a spectral type earlier than B9V. 


SU Gas: The binary nature of this Cepheid is still ambiguous. 
A component was first detected by Evans ( 1985| l based on study¬ 
ing the CAII H and K lines, and then a spectral type of B9.5V 
was determined from lUE spectra ( [Evans[[T991[ ). Although the 
location of SU Cas in a two color diagram is consistent with the 
presence of a companion, the radial velocity measurements do 
not show convincing evidence. [Szabadosj ( [199 l| l found four pos¬ 
sible orbital periods from the observations available at that time. 
Later, ( [Gorynya et al.|199^ compiled more data and suggested 
SU Cas as a possible spectroscopic binary. They derived an or¬ 
bital period of 408 days with an eccentricity of e = 0.43, but 
[Groenewegen[ ( [2008| l could not confirm the eccentricity with a 
larger data set and found a period of 407 days, fixing e = 0. Re¬ 
cently, Evans et ar| ( |2013] l reanalyzed special dates of Gorynya 
data, where the velocity difference was supposed to be the largest 
from the derived orbit, but they concluded that the orbital motion 
from radial velocity data do not show any convincing detection. 

Using the distance d = 392 pc (from the /T-band P-L relation 
of first overtone pulsators of Bono et ar][2002 1 and a Cepheid 
magnitude H = 4.27 mag ( jCutri et'H? 2003| l, a B9.5V com- 
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Table 5. 3cr average detection limits of the flux ratio. 



All Observables 
r < 50 mas 

r <25 mas r > 25 mas 

Only CP 
r < 50 mas 

r <25 mas | r > 25 mas 

Sp. Type 
upper limit 

Instrument 

VI334Cyg 
(2012 Jul. 27) 

1.71% 

1.45% 1 1.72% 

1.54% 

1.19% 1 1.57% 

B7V 

MIRC 

VI334Cyg 
(2012 Oct. 01) 

2.59 % 

2.15% 1 2.60% 

2.44 % 

1.86% 1 2.42% 

- 

- 

AX Cir 
(2013 Jul. 14) 

0.36 % 

0.40% 1 0.34% 

0.73 % 

0.80% 1 0.69% 

A5V 

PIONIER 

RT Aur 

0.64 % 

0.44% 1 0.64% 

0.47 % 

0.33% 1 0.47% 

FOV 

MIRC 

AW Per 

0.62 % 

0.54 % 1 0.62 % 

0.72 % 

0.52 % 1 0.73 % 

B9V 

MIRC 

SU Cas 

1.43% 

1.06% 1 1.44% 

1.37% 

0.99% 1 1.37% 

AOV 

MIRC 

T Vul 

(2012 Jul. 26) 

1.04% 

0.84% 1 1.05% 

1.13% 

0.81% 1 1.14% 

B9V 

MIRC 

T Vul 

(2012 Sep. 30) 

1.21% 

0.89% 1 1.21% 

1.29% 

0.92% 1 1.30% 

- 

- 


panion should give a flux ratio ~ 1.8% in the N band, which 
should be detectable with MIRC. Taking th e velo city amplitude 
of - 1 km s“' from Groenewegen (2008 while Gorynya et al. 
( 1996| l derived a velocity amplitude of 3kms '), we estimated 
a sin/ > 0.1 mas. Our nondetection could be explained if the 
companion was located at < 0.5 mas because it would not have 
been spatially resolved, but this type of close component would 
have an effect on the radial velocities (unless the orbit is face-on, 
but the probability of this kind of an orbit is low). The other pos¬ 
sibility would be that this companion has a wider orbit. We there¬ 
fore searched within a 100 mas radius range, but we did not find 
a significant detection. We then estimated the sensitivity limit 
for r < 100mas, and found a maximum 3cr flux ratio of 1.65 %. 
This means that if the B9.5V companion was within 100 mas, 
we would have detected it. Therefore, if this companion exists, 
it should have a wider orbit. 


The average sensitivity limits for a given radius range are 
tabulated in Table We reached a mean contrast of 1.37 %, 
which means that Aoth > 4.6 mag. This converts to an upper 
limit of the spectral type of an AOV star. 


T Vul; As with SU Cus, the radial velocities of this Cepheid do 
not show any signature of an orbiting companion, while a hot 
A0.8V component was detected by |Evans] ( | 1 992b| l . In the litera¬ 
ture, we found contradictory estimates of the orbital period from 
radial velocity measurements; for instance, Kovacs et al. (|1990| l 
found a long period modulation of 738d, [Szabados fl991[) es¬ 


timate d P - 1745 d from a larger data set, while BersieretaL 
( 1994| l did not find any orbital motion larger 0.55 km s“^ us- 
ing additional more accurate data. They showed that the long 
period of [Szabados (1991 1 is not compatible and argued that 


the 738 d period might be an artifact of the time sampling be¬ 
cause the observations were only made in autumn. [Kiss & Vink6| 
( |2000| l reached the same conclusion with additional measure¬ 
ments showing no signature of orbital motion in the radial ve¬ 
locity curve. 


From our interferometric observations, we did not detect any 
companion within 50 mas. The A0.8V component detected from 
lUE spectra would correspond to a flux ratio in // of ~ 0.7 %. 
According to our estimated interferometric detection limit, listed 


in Table|^ this kind of a component is below our detection level. 
We reached an average sensitivity limit of 1.04%, and we can 
therefore exclude any other possible companion with a spectral 
type earlier than a B9V star. 


5. Conclusion 

We presented an overview of CANDID, a new tool to search for 
point-source companions and estimate the sensitivity level from 
interferometric observations using the squared visibilities, clo¬ 
sure phases, and amplitude of the bispectrum, when available. 
CANDID allows us to; 

- efficiently detect companions using a grid of fit and deter¬ 
mine the detection level by giving the number of sigmas; 

- set the detection limit for a companion in data where a com¬ 
panion has not been detected; and 

- set the detection limit for a tertiary companion, in the case 
where a companion has been detected. 

We used CANDID to investigate a sample of binary Cepheids. 
We first determined the detection level for our previous detec¬ 
tions ( [Gallenne et al.|[2()13b| |2014| l and showed that the com¬ 
ponents were detected at > 13cr for AX Cir and > 25cr for 
VI334 Cyg. We also reported a new detection for AW Per, with 
a detection level > 15cr; the companion is located at p = 32 mas 
and PA = 67°, with a flux ratio of / = 1.22%. The compan¬ 
ion orbiting RT Aur might have been detected at 3.8cr, using 
only the closure phase signal, however, more observations are 
needed to confirm the presence of this component. Any addi¬ 
tional companions were not detected signffcantly (i.e., with a de¬ 
tection level > 3cr) around these stars. Likewise, no companions 
were detected around SU Cas and T Vul. From these interfer¬ 
ometric data, we were able to set upper limits for the spectral 
types; we found no components with a spectral type earlier than 
B7V, A5V, FOV, B9V, AOV, and B9V for VI334 Cyg, AX Cir, 
RT Aur, AW Per, SU Cas, and T Vul, respectively. 

The fitting procedure of CANDID also allowed us to mea¬ 
sure the uniform disk angular diameters of the new Cepheids 
observed. We found for RT Aur 0 ud = 0.699 + 0.011 mas (at the 
pulsation phase (j) - 0.32), for AW Per 6 ud = 0.627 + 0.018 mas 
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(at (p - 0.52), for SU Cas 0 ud = 0.609 + 0.043 mas (at <p - 0.77), 
and 0UD = 0.608 + 0.013 mas and 6 ud = 0.635 ± 0.018 mas for 
T Vul (at p = 0.27 and rj) = 0.12). 

We demonstrated that the approximation we used to analyti¬ 
cally inject a companion and estimate the detection limits is valid 
(i.e., error < 0.5 %) for contrasts /< 5 % if we use all of the ob¬ 
servables, and up to /< 50 % using only the squared visibilities 
and the closure phases. This makes CANDID a useful tool for an¬ 
alyzing long-baseline interferometric observations of binary star 
systems. 

Finally, this work demonstrates the capabilities of the MIRC 
and PIONIER instruments, which can reach a dynamic range 
of 1:200, depending on the angular distance of the companion 
and the (m, v) plane coverage. In the future, we plan to work on 
improving the sensitivity limits for realistic data through better 
handling of the correlations. 
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